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ABSTRACT 


Whenever active or passive devices are included in ar 
electronic circuit that is to be analyzed by & computer, 
appropriate models for these devices must be developed. 

The majority of the models that have been. adaptest foor ana-- 
dyzing a circuit on 2 cemputer have negtected tire effect: of: 
‘ambient temperature on the element. Temperature variation 
has been shown to have hazardous effects an electronic. 
components. This paper deals wrth the modificatiom of some 
existing models to include temperature effects anc stows: thes 
results of such modifications. It is done for six. bassc: 
PcG roir1e Commonents {1.€., BESIStOrT, Capacitor, imiductor, 
mode. bipolar transistor and field-effect’ transistor). 


Suggestions for extension of this work are also included. 
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i. UNTRODUCTION 


Ae GENERAL 

It 1s necessary for the electronic desigm engineer: to: 
predict and identify the response cf electromic systems: to: 
temperature effects caused by the device itself or: the: 
environment. Therefore, in the study of semiconductor: 
electronics the destrabirlity of computer simulatiom. becomes: 
readily apparent. The equations describing: semiconductor: 
behavior in all states are complex and nonlinear;: hences. 
solution by manual means ts difffcult even. wher many- 


Simplifying assumptions are made. 


B. BACKGROUND : 
1. Circuit Analysis by Computer 

WicCimemGempubter 15 USedyto andlyze am elecrronieG 
@mrcuit, a sufficient amount of information eoncerning the 
circuit parameters must be supplied and device terminal. 
Bewationships in the network equations describing the 
@ercuit must be determined. This is not a difficult task 
for passive components such as resistors and capacitors. 
Active devices may be described to computer analysis pro- 
grams in terms of models which mathematically simulate 
ethe behavior of the devices. It has become increasingly 
MmipOnaant to develope accurate madelLs of these devices: for 


computer programs used in design and analysis of such circuits. 


Analyzing the networks requires thet the derived 
equations must be in some particular form anc thear parame:- 
mers Specified properly by the user. The ciraurt elements: 
used in a network could include lumped resistances,. induc- 
tances, and capacitances, curremt and voltage sources: 
(dependent and independent), amd ective device models (i.e... 
diodes and transistors). 

In the application of the Computer Aided Aralysas: 
Peoprams the needvetten @rises for a circuit element that: 
is not included in the List of standard elements... The 
necessity of using linear but nonstandard elements: (such. 
as small-signal transistors or ideal transformers) arises: 
in many programs. In the case of linear nonstandard ele-- 
ments, general equivalent circuits have ofiten been developed. 
by others; they can be applied to programs. 

Zeer oimilation Of Cirecules with Active Fr llemenes 

ING Eas, SIS Kear SS; ele Nel\esliin lsiersueoleiie (bie Sicic slo! 
voltage sources or negative resistors. When we attempt to 
analyze an active circuit by computer we encounter some 
problems. The first one is to find a set of mathematical 
Equations to relate the current and voltage characteristics 
@recach device. These equations can be obtained by one of 
two ways; they can be derived from mathematical description 
-of the physical phenomena occurring during the operation | 
of the device or might be obtained by empirical observations: 


of the behavior of actual device. The CETTE oar ies 


ees -—= 





prediction, over a wide ramge of operatiom about artuad- 
behavior of the device, determines the value of these< 
equations. Secondly the parameters im these equations: must: 
be known or measurable or the required information. should. 
be derivable from terminal measurements.. Thee otheay factor: 


is the tolerance of the perameters.. 


Tl. ‘SBEVICE-MOMELEING METHODS 


A. INTRODUCTION 
The objective of modelling ts to proceed from the: 
fundamentals, that is, from matertal propertiess and mor-- 
phology of a distributed device te @ model as: presented: 
below: 
(1) To obtain systematically a lumped.,. een ee 
model for a distributed device, @) ame should writes 
equations which are distributed analogs: to ome of: 
Kirchhoff's two network laws, the current-node law 
or the voltage-loop law. Since the desired: model: 
is network-like, it must obey the two laws: of 
heerehhotf:s b) One should remove, by some method 
of approximation, all aspects of the mathematical 
Hescription which Can be traced to the distributed 
Mmeenre of the device. 
(2) To obtain models which describe the transient 
Penrormance of a device, ome might take advantage of, 
if possible, knowledge of its steady-state properties. 
The reason for this is, of course, that many aspects 
Omestéeady-Sstate beha@wior are much easier to déter- 
mine than their transient counterparts. 
In practice, sa medel is rarely perfect... Thus: we: 
Snouldsseareh fora model whose element= provide theabest- 


combination of linearity and frequency independence. 
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In this section the four basic approscies wroel be 
presented for developing equivalent cirawtts for elements: 
that are not standard in computer programs. These approaches: 
can be used individually or can be combined ta farm more 


complicated device model. 


B. PIECEWISE LINEAR MODELLING 

In piecewise linear modelling technique nominee 
elements can be analyzed by replacing all the nonlinear 
components with equivalent circuits, constructed with. 
mnie linear elements recognized by the computer program.. 
Tme nonlinear characteristics are averaged aver the rage: 
of interest and represented approximately by Linearized 
characteristics. For devices passing from one operating 
region to another (such as a transistor Pe aat from the 


Bee cerecTon £O Saturation), a different. linearized 


eivalent Circuit can be proposed for each recion. 


C. MATHEMATICAL MODELLING 

Those circuit elements that are not standard (defined 
mum Circuit element that 1s not include@€d in the Ent of 
built-in models of any computer library) can be modelled 
using only the mathematical relationship which describes 
miemin terminal characteristics. The procedure to follow 
in developing equivalent circuits using mathemattcal model-- 
ling technique is quite straightforward. Im many cases. 


the modelling problem can be reduced to the synthesis of 


iat 


a simple RLC network sometimes requiring tire usxe of 
dependent current sources. Im other ceses,, = little 
inventiveness is useful. 

The equations describing the components must be linear 
Since only linear elements are used in the development: of: 
the model. The availability of R, L, are C elements: im. 
the construction of the equivalent crreurt permits: the: 
modelling of linear differential equattans with cons:tant-: 


coefficients. 


D. PHYSICAL ANALOGY: PEECEWISE LINEAR MODELLING 

The modelling of a nonlinear element with an algebraic: 
voltage-current relationship (such as nonlinear resa's:tance:) 
mas discussed in Section B. dm this section the techniques 
for modelling a nonlinear reactive element with a non- 
algebraic voltage-current relationship is discussed (such 
aoeinauctance or capacitance). For better understanding 
Eiicmmay reter to €xamples given in Ref. 1/7. Since non- 
linear inductance and capacitance are algebraic functions 
Giemacnetic flux and charge, respectively, which are not 
recognized by the program, a reasonable solution to this 
modelling problem is to find a physical analogy which will 
rglow flux or charge to be treated in those units: recog- 
mezed by the program. 

The physical analogy acts as a black box which converts 
ene physical reiatmenship into terms witich are acceptapie. 


Eom tne program as Shown in tollowing algorithm. Therefore, 


Fe 





elements such as nonlinear imductamces and Capecttarces: cam. 






PHYSICAL PHYSICAL COMPONENT EQUIVALENT ss 
CIRCUIT USING COMPUTER 
COMPONENT ANALOGY V-I RELATIONSHIP PROGRAM 


Figure 2-1. Physical Analogy Cancept Use 
in Physical Component Modelling... 


be modelled if an analogy between the characteristics: to: be: 
modelled and a voltage-current relationship acceptable to 


program can be established. 


E. PHYSICAL ANALOGY: MATHEMATICAL MODELLING 

Almost all of the Computer Aided Circutt Analysis 
Programs perform calculations in terms of voltages and 
currents only. The parameters, such as charge and magnetic 
mite. have not been included in the programs and Cannot De 
mer-amatrectiy. if Such Conditions Exist, 1t 1s mecessary 
Gemdevelop methods for expressing the new parameters (charge, 
meemetic flux, etc.) as functions of voltage and current or 
to redefine some new elements in addition to the elements 
included in the standard list that will perform the required 
analysis in terms of voltage and current. 

The solution of mechanical and thermal problems by 
electrical analogy techniques is an example of this: kind. 
Establishing the analogy is the first step in the deriva- 
feon Of an equivalent circtmt. The second step 1s: linking: 
the derived analog network to the remainder of the model. 


which includes the element terminals. The element terminal 


53 





characteristic must always be defined im tenms of voalltage 
and current. The charge-control tramsistor model ts & 
good example of physical analogy techniques tn modelling 
the other charge equations which govern the tranststor 


mesponse. See Ref. 17, Appendix A.Z. 


Jeu 





Mri TEMPERATURE EFFECTS 


fee DEMPERATURE EFFECTS ON RESISTOR 

The resistivity p of a conductor is the ratio of the 
Peet ric intensity E to the current per unit cross-sectional 
aviee a. 


p = E/(I/A) Gail 


Mnremresistivity of all conducting material is affected by 
Mmemmcenperature. A plot of resiStiVity versus temperature, 
mmeinctallic conductor, 1S given in Fig: 3.1. The curve 
may be satisfactorily represented by an equation of the form 


2 3 
a Geet oy Be ee eae, : ome 


where Py Hone wresistiy Ley eateul: 1. 84> g>, nat, eee 


constants depend on the characteristic of a particular material. 


—_ 


a Se 


boure S-lo me ReSIStivity versus Temperature 
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For temperatures which are mot too great the termss im 


2 


T” and higher powers may be neglected and cam be written ass 


Pp = Op + git : be: 
The temperature coefficient definition of resistivity: is: 
a= g,/9 = (o-0,}/Co,T) 35, A 

Since the resistance of @ given comductor ft proportiorad. 
to its resistivity the following equation can be written 
R = Ry + saa ) Dice: 


where Ro is the resistance at 0° C and R fs the resasstance 
oer =«(C. 

The resistivities of nonmetals decrease with increasing 
Memperature, and their a's eresmerative. Electrotytes- also. 


4 


have negative a's. 


PoeeeMPERATURE EFFECTS ON INDUCTOR 

Temperature variation is usually accompanied by a change 
Pueemauctance, resistance and self-capacitance of a coil. 
tie last two types of change are not important in. compari- 
weeewith the change in inductance. The formula for sing ile- 


jJayer solenoid given by Sturley [Ref. 6] is 


L = (r“N*)/(9r + 102) 3.6 


where & = length of winding 


be) 
i 


radius of winding 


z 
I 


tetal turns  iimene "cor! . 
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The change of the inductance fom am increase im temperature: 


of T° C can be skown to be [Ref.. €] 


=" b GA) oa 
where AL= change in inductance 
b = coefficient af Linear expanssom of thee conductor: 
AT = change Im temperature.. 
Equation 3.7 is valid whem the coefficremntss of: radiad: and: 
axial expansion, b and D.,.,. are approximately equal.. 
Htets possible to make the inductance. L independent: of: 
eempecrature by satisfying Eq.. %..&,. the: ratio: of: axiad: toc 


radial expansion. 


b,/b_ = toate, il) 2) ory oss 


By choosing the coitl-former and the conductor material 
jMeeperiy, the temperature coefficient for inductors can be 


minimized, 


Wee TEMPERATURE EFPECTS ON CAPACITOR 

Temperature variation has two effects on capacitors, 
namely to change the dielectric permittivity and the dispo- 
wrenon Of the plates. sSeme typical capacitanee-tCemperature 
coefficients will be given in Chapter IV. Permittivity of 
the insulating material is considerably changed with tempera- 
mine. Capacitances for parailel plate-and cylindrical types 


are given with the following equations, 


parallel plate: C= (A.ej)/d 339: 


ey 





cylindrical: C = (2m. S.e)/lTog (a/b) 38 100 


where area of plate (A), and length cf capacitor (‘Q)),, are 
temperature dependent also, but this dependence can be 
negligible when compared with the dependence of permittivity. 
In both cases capacitances are the functtonm of permittivity 
ce, In addition to temperature dependence af ¢,. it iss 
dependent of frequency, but the Latter effect is a different: 
SPOry . 

Improvement of the manufacturing techniques mimimizes- 


the temperature coefficient far capacitors.. 


fee TEMPERATURE EFFECTS ON DIODE 
ite introduction 
Mnere sare seve Kiidsiot GGitaigie —G amiayal ihe: scat taleciic Sa 
one 1S negative (the free electron) of mobility HL» and 
meme other is positive (the hole) of mobility Uns These 
particles move in opposite directions in an electric field 
Pemout Since they are of @pposite sign, the currents are in 


the same direction. Hence the current density J is given by 
J = (nu, + Pu) e E 


where n = magnitude of free-electron concentration 


Pp magnitude of hole concentration 


e electron charge in coulombs. 
For an intrinsic semiconductor n = p = Ns». where n, is: the 
mitrinsic Carrier Concentration. he intrinsic carmcer 


concentration n. 1s temperature dependent, such as 
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Zoe 3 : 
iA AT~ exp ( Ego/ kT) 


It has been shown experimentally that Eo depends upon 


temperature also as pointed out below 
Eo = Eco = gi 


where Eco is the energy gap of the semiconductor at 0° K 


4 


and q is a constant (q = 2.4 x 10 ° eV/°K for silicon and 


q = 3.9 x 10°" 


eV/°K for germanium).* Such a large change 
conductivity with temperature places a limitation upon 
the use of semiconductor devices in some circuits. 

2. Temperature Dependence of Parameters 


For a p-n junction the current I is related to the 


voltage V by the equation 
I = I, lexp(V/mV.)-1] Sel 


iemperature dependent parameters of the Eq. 3.11 are Te 
Voltage 0 








1 Micro amp. 
Bo BeS(C 

LOW eiite ro ain ae 
COMBS oie 


(b ) 


Figure 3-2. Typical Reverse Biased Diode (a), 
anid mes Cm@ancetenr is pres (0) 


motcrictiy speaking, q_is not a constant. For example, for 


seirecan q = 2.4 x 0 eV/°K only when T > 200° K. For 
germanium, q = 3.9 x ie eV Komik whem, dS 2 ic, 


ES 





i is called the reverse saturation qurremt ami eivear by’ 


the equation 


= n = - a 
V3 = K T exp( Veo/mV 5) Sr bes ies 
for Germanium: m= 1 m= 2 Veg = 0..7&5 valt: 
mor Silicon: m= 2 m= 1.5 Veg = ..20 volts 


From experimental data we find that the revemse saturation 
current increases approximately seven per cent per * C.. 
Sence gine? )ae = 2.0 reverse saturation current approximate:- 
ly doubles for every 10° C rise in temperatune.. Va Iss thes 


“wolt equivalent of temperature” defined by 


Ve = T/11600 . 


mmotier important parameter for dai@de 1s diffusiom ors-torage: 


é 


Sageacitance and given by 


Go = Galina cle 


D T 


where t is the diode time constant equals the mean lifetime 
Pamiinority Carriers. 
meverse leakage resistance also Nas some temperature depen 
demee; Since it can change depending on the manufacturing pro- 
mess, nO generalized temperature coefficient can be stated. 
It should be determined for each family of the diode. 
maectempenature coefficient of the diode parameters can 
bemacrived in differential form from Eqa? 3.11 and 3.11.4 
by taking the derivatives and making some approximations 


toryield the following equations ¢ 
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(1/I,} (41 /dT) = (a/T) + (Ugp/eTV,) 35 15 
dV/dT = [V-(Vgg * mV) ]/T G8 bee 
at fixed current. 


E. TEMPERATURE EFFECTS ON TRANSISTOR 
1. Introduction 
A junction transistor consists af a salicom or: 

germanium crystal, in which a Layer of n-type or p-type: 
material sandwiched between two layers of p-type or: n-type: 
material, is called @ p-n-p or n-p-—m tmamsustorm. Fort as 
@ualitative umderstanding of the form of the imput: and: out-- 
put characteristics, ome cam consrtder that the transistor: 
eemei1sts of two diodes placed rn series" "back to~baek." “It 
has four operating regions defined as follows: 

meerom (a)28 CUuL-Off region (emitter reverse biased, 

collector reverse biased), 

Region (b): Normal transistor operation (emitter 

forward biased, collector reverse biased), 

maerome(c); Scaturatronm=negion (emitter forward biased, 

collector forward biased), 

Region (d)}: Inverse transistor operation (emitter 

reverse biased, collector forward biased). The 

mIpuc-OUtTDUt Characteristics of the transistor in 


different modes are given im many texts... 


call 





2. Temperature Dependence of Parameters: 

All of the parameters af iInteresst br computerr 
application concerning transistors change with temperature. 
Some of them are very sensitive to temperature. changes., 
while others may be assumed essentially coms:tamt: im the: 
analysis. The most temperature sensirtive parameters; ares 
Saturation currents, Tec and Lies: current gaanr. factorss. 

Oy) and Ars and leakage resistances,,. Ry and Roe. Parameters: 
like time constants, Tue Type capacitances C..,. C.,. propor>- 
tionalities, Mes Mes change only slightly. over: the: 
temperature region of interest as high ass 125°° C. 

(a) Reverse saturation currents: The emitter: 
junction diode saturation current [zg and collector: junction 


diode saturation current I are the mMOSt .bempemature Sensi 


cS 
tive parameters. Suppose that the reverse saturation 
current at room temperature (300° K) is given by Eq. 3.11l.a; 


it can be written at temperature T + AT, as follows: 


len = 6) T exp l-EeaeaT) 
ee eNe peptcEen/ s(t # 277) 3.15 
equivalently 
Pao. = Katt + AC Ne exp[-Eco/k(T + AT)] 3.16 


Esoumang AT << T and expanding the Eq.. 3.15 in Taylor series, 


the following can be obtained: 


Zz 





Teor = Icol) *iCEgg/RF) * mE (AT/T)P = Leg lil + MCAT) 3317" 
After numerical substitution it turms out to be 


I = Teg GQ + MO iog GAT) 


COT 


where i» = 0.157/°C and is defined as the temperature 
coefficient for silicon tramsistor. For greater accuracy). 
over larger temperature ranges, the equatiom tncluding: 
mre Exponential factor directly should be used. 

(b) Current gain factors: The current gain. factors: 
in normal and inverse modes for common base configuration 
are given below: 


Qy = hpp/(1 + App) and a; = Beey/ (1 + Degy) 


where Der is normal mode, common emitter current gain,. and 
Ape] is inverse mode, common emitter current gain. oN and 
a do not change appreciably with temperature but (1-a,,) Om 


(l-a,) plays an important role in the following formulas: 
Dep = dy / (1-a,)) and Dery = ar/(1-a,) ; 

Fer example, the variation of common-emitter current gain 

hops wrth temperature normalized to unity at T - 25°C 

memeche 2NSZ4 transistor 1s shown in Fig. 3.3. 

(c) Reverse leakage resistance: Reverse leakage 
resistance is temperature dependent to some extent. No 
memecralized temperature coefficient can be stated sinte it 
Vortocmudenendine Won thegpmacessuat manutacturings, There-- 


fore, it should be determined for each family of transistor. 
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= 0 - 60 -20 20 60 100 
Temperature °C 
Figure 3-3. Dep versus Temperature 


F. TEMPERATURE EFFECTS ON FIELD-EFFECT TRANSISTOR 
fie introduction 

Preld-effect transistor (FET) has two primary 
operating regions, namely the triode region in which the 
device behaves as a voltage-variable resistor, and the 
pinch-off region where the current is determined mainly by 
the gate voltage. The two regions are separated by the 
ies Of Pinch-off points, 1.e€., the points where Vos = 


Voee- =v V ea iG 


D eon DSee P 


menecace, pinch-off voltage and gate to source voltage, 


Gs are defined drain to source 


meopectively. 
the operation of FET depends upon the flow of majority 


marriers only. Therefore, it 1s a unipolar device, since 
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the majority-carrier current decreases with temperatures. 
the troublesome phenomenon [Ref. 4] Is not enmoounterecd: 
math field-effect transistors. Whis property provides: 
thermal stability for FET. 
2. Temperature Dependemce af Parameters 

Drain Wha fron Cteereiit,. Ines Which depends: upon 
channel geometry, impurity density and charge carrier 
mobility, 1s defined to be the channel current,,. Dys3: 
at Ves = “Vp. Although Ingg Is mot greatly dependent. on 
Vos: it depends on temperature heavily... Temperature cam. 


influence the drain current In two ways.. Fimstly,, a changes 


in mobility causes a change in resrstance of channel... 


‘Secondly, a change in temperature causes: a change im. surface: 


Mercntial. The chamge in surface potential warks: lr opposa-- 
tion to the change in mobility. Thus, change in surface 
Metentcial constitutes a built-in gate-voltage compensation 
faeene mobility change in the semiconductor. The effect is 
dependent on the D.C. gate bias so that at a particular bias 
the compensation may be optimized. 

Another quantity of interest rs the transconductance. 
mie) transconductance, at has the same temperature variation as 
Ine does. The principal reason for the negative temperature 


egerricient of I is that the mobility decreases with 


DS 
mrercasing temperature. 
AS POM teCdmOnE Deteone, (Eine) Magon Citecits- on tempemda. 


mire On a FET or aiteration of the average channed conductivity 


LS 
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DSS 
es, 


= 50 0 50 100 SO. 7 
Figure 3-4, Normalized Ines versus Tempera- 
bbe ee se 


mmecne gate-channel junction barrier potential. Since 
memductivity decreases with increasing temperature, the 
drain current decreases as the Ce eee On the 
@ener hand, barrier potential decreases with increasing 
temperature, causing the magnitude of the total gate-source 
Peeeseral to be reduced, hence increasing the drain current. 
By suitably choosing the operating point, it is possible to 
make these effects canceled, producing what is essentially 


mmeeno temperature coefficient operating point. 
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IV. COMPUTER MODELS 


A. INTRODUCTION 

Almost all the computer analysis programs are destgned 
to solve the linear and nonlinear simultaneous equations 
which characterize the mathematical models used: to predict 
circuit response for an electrical/temperature environment... 
iimemser codes the problem by entering the circurt topolagy 
and parameters. Programs contain 42 standard (built-in) 
component model library. Non-standard models may be can- 
Structed by using standard models. FORTRAN algebraic 
equations and control statements may be coded to modify 
standard model parameters, augment the standard models, 
create nonstandard models, create additional time functions, 
Poeemature functions, etc. The cirevit 15 entered into 
the computer program by using the program standard part 
models and specifying the circuit topology. 

In this paper by "conventional (or non-standard) model," 
it means the model of a device which is isolated from all 
environmental effects, and by "modified model," 1t means a 
model with temperature effects taken into account. In 
addition to this task, some additions or some simplifications 
will be done when necessary, depending upon device characteris- 
ics . 

In order that electronic components may be handled 


quantitatively, a set of mathematical relationships must be 


Za, 





derived relating their terminal characteristics. Deriva- 
‘tions will not be attempted in this paper since they are 
documented elsewhere [Refs. 2, 3, 4, 7, 18, 22]; only the 


results will be presented. 


Pee RE oLSTOR MODEL 
1. Conventional Model 
The resistor model is defined by the following 


current equation 


pk —_—— Ao T 
1 R 2 


Figure 4-1. Resistor Element 


Since the element is bilateral, the choice of which is con- 
sidered to be terminal j] and which is terminal 2 depends on 
the direction of current flow (from T) to T.) . This element 
Nas juSt one parameter: the resistance value R. 
yee Modified Model 
Temperature dependence of resistor was discussed in 
mmaoeer Jil, but the Eq. 3.5 can be réwritten in slightly 


different form such as ; 
R. = R £+ AR Aa 


where AR is the amount of resistance change due to tempera- 


ture change and equivalent to 
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AR = +aR(AT) Artz 


The sign depends on material properties. retnine esas Hon 
aS linear, the following modified resistor model for 
Gomputer can be used. The determination of AR requires the 
parameter a, which can be found in many manuals or reference 


books. 


(IR oe mm Na NN) 


R AR 
Figure 4-2. Modified Resistor Model 


C. CAPACITOR MODEL 
ieee COnventional Model 
A capacitor, shown in Fig. 4.3, is a bilateral 
element with terminals Ty and T,; havine=a\ Capac teanecmae- 


The current flowing through a capacitor is given by 
1, = C{d(V, - V,)/dt] 


See | 
C 


ce | are 


Pigure 4-3. Capacitor Element 


Zo 





2. Modified Model 
The more realistic capacitor has a shunt leakage 
resistance and a series resistance as shown below that are 


ma@emtified with Ry and Ro; respectively. 





Figure 4-4. Capacitor Model with Resistive Components 


Temperature dependence of capacitor was discussed in 
C@apter III. Very small series resistance (in the order of 
0.1 ohms) and very large shunt resistance (in the order of 


a 


10° ohms) can be assumed constant with temperature change. 


some typical capacitance-temperature coefficients are given 
in the following table and many others can be found in the 


references given in the footnote.* 


wae 4-1,  Capacitance-lemperatune Coetirtecrents 


Insulating Material Cabal ane clei e tect disc mee etm neene mines 
Varnished Cambric +2100 parts in 10° per degree Centigrade 
Synthetic Resin +1600 parts in 10° per degree Centigrade 
Enamel +470 parts in 10° per degree Centigrade 
Ceramic +100 parts in 10° per degree Centigrade 


*(1) Handbook of Chemistry and Physics; (2) International 
Serercal fables Vol, 2, McGraw-Hill; (3) Tables of Dielectric 
Materials, Vols. I-IV, MIT, January, 1955. 
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Defining a change in capacitance as AC the modified capacitor 
model can be obtained shown in Fig. 4-5, in which AC should 
be calculated and added to the model with the aid of tables 
(ven in ie mmetnote on page 30) and the Eqs. 3.9 and 3.10. 


The change in capacitance is very important if the capacitor 


meoure 4-5, Capacitor Model with Temperature Effect 


is an element of resonance circuit because a small change in 
4 


capacitance may cause a large drift in resonance frequency. 


D. INDUCTOR MODEL 
1. Conventional Model 


mie relationship for Current 1s given by 
1, = (1/L) fs (V, = V,) dt 


where L is the value of inductance in henries and the dif- 
ference, (V, - V5), is the voltage across the element. The 
inductor is diagrammed in Fig. 4-6. Since the inductor is 
bilateral, the choice of which to be considered T, 1s com- 
pretcecly arbitrary. This element has two parameters: induc- 
tance L, and Q (of the coil). Many applications of inductors 


have practical equivalent circuit that also includes the 
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Figure 4-6. Inductor Element 


Memeity QO--not just pure inductance. Frequency dependent 


mers civen by 
i SO 2yaell A 


where f is the frequency at which the analysis is being 
mamcomlced, Ihe significance of Eq. 4.3 1s best seen in a 
Sketch of resistance versus frequency. The range for the 


Beemieael 21S that r becomes unrealistically small as frequency 


4 


r 
Band of <9 g 
reasonable K— 
;modelling 
| 
, : US oo ty DaGavleby. 
} specified at 
some test 
frequency. 
| 
| 
. fe 
0 


muotre 4-7, Resistance r versus Frequency £ 


becomes very low, and at high frequency the model would show 


meapproaching infinity which is also unrealistic. 
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Z. Modified Model 
Temperature dependence of inductance was given in 
@mapter Ill with Eq. 3.7. Additionally D.C. series resis- 
tance and the shunt leakage resistance should be involved 
into the model, which are assumed to be independent of 
temperature. The schematic which takes all the effects 


explained above into consideration is shown in Fig. 4.8. 





mieure 4-8. Inductor Model with Temperature Effect 


To determine AL one should first determine the parameter 
b in Eq. 3.7. The value of b can be found in the handbooks 
feeaed if the footnote in prévious section. For better 
mrcerstanding of importance AL, look at the example given 
below. 

Example: For a tank circuit, f = 1/(2nvLC). If Af is 
the Change in oscillation frequency due to AL change in Ly 
the following relationship can be obtained after proper sub- 
Srerution and expansion by the Binomial Theorem where only 


the linear term is retained: 


= re 


5p 





The value of 6b for capper ts L..6 x TO” pets: pear °° CC. Thus: 
the ratio change in frequency due to the unhindered: expansion= 
G 


Sia copper Goil is # = IG ~ pawts per ~ C.. or for: AT == 39°'C: 


the ratio change of frequency is 240 parts in 10°,. Frequency 
variations due to cotl expanstom of 240 amd 2400. hzz ati oseil-- 
lation frequencies of I amd 10 Mhz,. respectively.,. womldcbe<- 


expected under these condittions.. 


E. DIODE MODEL 

-1. General 

Any diode can be modelled. using the following: approach. . 

The quantities of interest are electrical charatteristicss 
which can be observed at the external terminals: (lLes.,. a 
Meerack-box'' approach). This type cf modeling requires: 
finding a Mathematical expression to errucae Sele Elec igr i = 
@aieecMaracteristic. In this paper this method will. be 
emphasized. Physical constants will be used whenever possible, 
but if a characteristic cannot be related depsiest ae an 
emparical description will be used. 

Based on the operating range, diode models may be 
classified into two groups: "small-signal" and "large-signal" 
models. The first type consists of finding an operating 
point and assuming that all responses are linear about that 
point. Therefore, this model 1s valid over only a limited 
operating range, mamely, where the response characteristic 
is linear. Large-signal model its one which anv eee 


a wide range of operation. By its nature, it is the best: 
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suited model for use in general computer analysrs.. Theres 
are three most widely known models: (1) EbersMalll [Ref&. 1D],. 
(2) Beaufoy-Sparkes charge-control model [Ref.. IV},. and 
my) Linvill lumped model [Ref. 12]. In this paper Ebers— 
Moll model will be examined im more detarl and necessary 
modifications will be given. The reason for this preference: 
is this model describes the "black-box" behavior with reda-- 
tively simple mathematical equations. 

Ze eModel withowt Temperature Effects 

The current equation for dtiade was gtvem before 1m 
Eq. 3.11. Another diode curremt equation will be obtained: 
by discussing the various factors in the semtconductor diode: 
which have to do with its electrical performance... 

> Ven the applleGcattoOnm Of al ex tcinal Tevense— ras 
eye p-n junction width will expand because more donors and 
meeemtors Will be uncovered, Ihis increase in uncovered 
charge with applied voltage may be considered a capacitive 
effect, named transition-region capacitance. 

If bias in the forward direction, the potential 
barrier at the junction is lowered and holes from the p-side 
enter into n-side. Similarly, electrons from the n-side move 
mimes p-sede. A change of voltage on the p-n junction will 
mbcer the amount of stored carriers. Thrs, then looks like 
aeeapacitor to the accompanying Circuit, named storage or 
diffusion capacitance. 

Byetaking the Leakacle Tesistance ito, account, the follow] 


mroecdiode current Equation can be written as 


HS 
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I, = I + E/R, + Cj, (dE/dt) 1" C..(dE/dt) 4.4 


These previously derived properties are now all combined to 
provide a composite equivalent circuit, or computer model, 


Geeethe p-n junction diode, shown in Fig. 4.9. 


Anode Ry 2 Cathode 


Figure 4-9. Diode Computer Model 


tae definition and determination of the model parameters 
will be discussed o the following section. 
3. Determination of Model Parameters 
Ideal model parameters can be read directly or 
extracted from the device Specification sheet. Unfortunately, 
this cannot always be done. Most data sheets do not contain 


sufficient information to define the entire model. 
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Bie ee Rp represents the semiconductor bulk ressis:-- 
tance, Ry the junction leakage resistance, C,. the transittomr 
Capacitance, and I is an ideal current source whitch describes: 
mire fLOrward current as a function of junction voltage E.. It. 
1s possible to write the Eq. 3.11 in slightly different form 


as given below: 
I = I [exp (eE/mkT)-1] 1 hs 


where e, k, and T are electronic charge, Boltzmanm's constant,. 
and absolute junction temperature in “°K, respectively.. I. 
and m are constants which define the shape of the V-I clrarac-- 


teristics. For D.C. operation, the network equations are:: 


Vy a * In Rp ¥ 6 
Ty = J se jo, 
A plot Of the static Vy - ry characteristic for 


positive In ie eSlvoyipe ai Ieee ole File lhe 1S leu ues 

on a logarithmic scale to illustrate that over a wide range 
mene forward operating region, Vy = E. For determining 

the unknown model parameters Ij, m, and Rp one method is 

to curve-fit data points on the forward diode characterrstic 
momeads., 4.6 and 4.7, where In 1s assumed to equal I. Another 
approach to finding the parameters involves the same assump-— 
tion but not the curve fitting. One may choose two points 

on the linear portion of the curve in Fig. 4.10 and label 


them (I Vi) and (I,, V5). In this range may be expressed 


1 led 


oy 
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with following: 


where 6 = e/kT. 


=, 1;.Y,) 


~“ 


Figure 4-10. 


Lac") 
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Diode Forward V - 


exp (6V,/m) 
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Taking the logarithm of both sides and substituting 


(I; Vi); (I,, V5); then solving for m gives 
m = 8(V, - Vj) /in(I,/1,) 2 4.9 


Substituting either data point and the value of m found in 


bape 4.9 into Eq. 4.8, the following can be obtained: 


I 
O 


I, /exp(8V,/m) = T,/exp(6V,/m) 


The difference between the linear extrapolation and the 
Measured point is the drop due to Rp: solution of Eq. 4.6 


results in 


Ry = (Ne SET 4.10 
where Ex = (m/6) ln [(1. + ted = (m/6) In (I,/1.). 
In 
-10 -5 0 
| 
! Z--+ to 
(1,,V,) -2uA 
Figure 4-11. Diode Reverse SOs 
V - I Characteristic 

-4nA 


From the reverse characteristic shown in Fig. 4-11, the 


Ohmic leakage resistance, R, can be found to be 


So 





Ry = V,/ Cy, + I.) : ie el 


There remain two more parameters to be determined, namely, 
the junction transition capacitance Co and diffusion capaci- 


tance Cy IMP oN ey ptcal Curve Of JUNCTION Capacitance 


Wersus reverse bias is shown in Fig. 4-12. 





inne 4-125 Junmetion l[ransi tion 
Capacitance Characters tic 


the analytical vexpressionm that deseripes Ci 1s 


Cae! R/O” “4.12 
where K, n, and V; ane JUNCCION Capacitance proportionality 
Gemstant, grading constant and junction contact potential, 
meepectively (V5 is approximately 0.9 or 1.0 volts). For 


V, = 0.0, C,= C= K. Therefore, at any point (C), vas 


— as 
Cc, SEU V1) : 4.13 
wemevne Eq. 4.13 for n results in 


n = In(C,/C,)/1n(V, - V4) 
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Another way of expressing the Eq. 3.12 may be the following:: 


Cy = @(I + I) smart 44, 15> 


where the parameter f is diode cut-off frequency and evalua- 
tion of f can be done by Kuno's* method. The cut-aff 
frequency is given by 
f = l/2Znt . 
4. Modified Model 
The temperature dependence of diode was: examined im 
Chapter III. Temperature coefficient of saturatiom current 
and temperature dependence of the diode voltage were giver 
Saoeg. 3.13 and 3.14. Equation 4.15 is e#pressed in more 
convenient form as 
Cy = O17 me Cewy iie ale aly) ee cane (46s) 
Pemesimplification purposes, making some definitions will be 


useful at this point such as follows: 


A 


I 


et/mk Ae IL 
B 


eE/mk 4.18 
G = Fractional imereds cme esGe ies 
current with temperature. 
Temperature dependent reverse Saturation current increasies 
exponentially as indicated by the following equation: 
I = JL ee(= 9.40) ROG e += La 4.19. 


oT O 


fae J. Kuno, “Analysis and Characterization of p-n Junction 
intede owitching,’’ [EEE Trans. on Electron Devices, Vol., ED-11, 
pe o-14, Jan. 1964. 
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where lor 1s the reverse current at temperature T: amd [iss 
the reverse current at temperature To The solutiom of 


Eq. 4.19 for C turns out to be 
C = In(I ./1,)/(7 - TJ 2 4-, 2D. 


For example 1N3605 silicon diode has the following: specafica-- 


tions: 


—_ = 1°] 
I, =0.05 mA at Ty = 25° 9G 


Applied Voltage = --30 volts. 


— = 1°] 
lor = 50 mA at T = 150° C 


Therefore 


C = 1n(50/0.05)/(15@ - ZS) = @..055/ °C 


Proceeding with concise form of Eq. 4.16, the following 


relationship is obtained: ‘ 


Cy) = A(I/T) . 4.21 


If we drop the unity in the parenthesis in Eq. 4.5 an approxi- 


mate current equation can be written as 
I = I fexp(eE/mkT) |} = Is exp (B/T) beg ZZ 


where B is given in Eq. 4.18. The AT change in temperature 


causes the following modification for Eq. 4.22. 


le 


7 = jp exp(B/(T + AT)] 


ie exp(C.AT) exp[B/(T + AT)] 


Ie exp(C.AT) exp(B/T) exp[-B(AT)/T(T. + AT) ] 
Definition of another factor such as 


D = exp(C.AT) exp[-B(AT)/T(T + AT)] 
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provides a very short expression for I given below: 


T? 


I. = Dae EAS 


On the other hand AT change in temperature implies the fol- 


lowing modification on diffusion capacitance C 


i 
Cop = Bolle = BID) Ss NNR Ie Gr ay 

= (A.I/T)[T.D/(T + AT)]= Cp) . T.D/(T + AT) 

= F.C 4.24 


D 
pierce F = (T.D)/(T + AT). Then adding and substructing the 


Cy) to Eq. 4.24 yields 


C Ba CAG Ca 


DT Sip) D D 
= Cy + Cy (ie le), Aes 
In similar fashion 
I. =. Del. 7] 
Sale ce wi Gb) 4.26 


The temperature treatment of more sensitive diode parameters 
was given by Eqs. 4.25 and 4.26. The changes due to tempera- 
ture for Rp and Rp 


Applying all necessary modifications one would get the com- 


are assumed very small and neglected. 


puter diode model when temperature vatiations exist as shown 


eig., 4-13. 
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I @eo) i 


Cathode 


Figure 4-13. Diode Computer Model with Temperature Effects. 


Hee RANSITSTOR MODEL 

meee oeneral 

A transistor model in computer-aided analysis should 

meet the following properties: (a) The model must consist 
emeclrcuit elements that are recognizable by the computer 
program being used. (b) The parameters used in the model 
Smemra either be obtainable from published data or readily 
estimated. (c) The model should be valid to high frequencies 
as well as to low frequencies. (d) The model should function 
without regard to the circuit configuration (common base, 
common collector, or common emitter) in which the transistor 


moeplaced. 
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Note: A general model for a device does mot always lead toc 
the most efficient analysis simce certaim troublesome elements: 
of the model may not actually be required by the problem at. 
hand. For example, nonlinear or high frequency elements: 

may not affect the circuit response greatly, but could 

easily affect the solution time significantly... 

There are many models that have beem praposed to: 
represent transistors at low frequencies. Any af these small-- 
sremai models 1s valid only when the transistor operates: 
in the active region and produces negligible distortion 
in the output signal. Fortunately the output characteristics: 
Dietransistors are reasonably Linear aver @ fairly large 
range, often allowing output swings of several volts with 
Mictle distortion. High-frequency models for the transistors 
are more complex than low frequency models, since additional 
Metive elements are included. 

eee Model without temperature Erftect 

The transistor is a three-terminal device, fabricated 
Pyeolacing two p-n junctions in close proximity to each other. 
The current through one forward-biased junction strongly 
imordences the current through the other. The normal mode 
Meeeransistor operation occurs when the base-emitter junction 
is a forward biased and base-collector junction is reverse 
biased. <A transistor may also be operated in the inverse 
mode. In this case, the collector is forward biased and: the: 
emitter is reverse biased. Therefore, transistor can be. 


considered as two junction diodes connected in Series, of. 
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which one is forward and the ether is reverse biased... The 
emitter junction and the collector junction current equa-- 
tions can be written as in single diode case. Sy referring: 
to Section E and taking the inverse mode current gain into. 
account, the following current equations may be expressed in 


terms of diode parameters: 


I; Ter i (VE/Ry gE) + Cop (dV. /dt) + Cy (dV / dt) = oc -L 4,27" 


In = Top oF (VO/R, ©) Cro (dV ./dt) is Cy (dV./dt) - OT 4.28 
where parameters Oa: and a, are mormal and tnverted short-- 
circuit current gains between emitter and calltector,, capaci-- 


tors CoE and C fp represent emitter-base depletiton—lUayvexr aad: 


D 
emitter-base diffusion capacitances, respectively... White C.u- 
and Cc wepresent the collecton— pare aden le tron cmc edsiel tl oamom 
capacitances, Ve and Ve ane tne emitter base pana weelle Gear. 
mee vVOltages, respectively. Emitter and collector reverse 
Saeetraction currents may be involved into the junction current 


equations in such a manner given in the following equations: 


Tap = I.colexp (6V_/m,) - I] 4°, 29 
Top = Tog lexp (8V./m_) = Wale 4 ..30 
Tic and Lac are usually given in specification sheets, but 


the analytical expressions for Tes and ers are given as 


ES 
and 


cs = Fag/ ~ OO) 4°,.32- 
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where I is the saturation current of emitter jumcttiaom wit: 


EO 
@ero collector current and Tao is the saturetion current af 
collector junction with zero emitter current. 

Now we are in a position to model the transistor 
device. The computer model uses two p-n junctions modelled 
like the diode previously discussed plus the current ga@ir 
factors as shown in Fig. 4-14. The resistances Re and R~ 
account for ohmic effects and can be approximated as constants 


in most instances. 


‘ The junction capacitance factors are: 


: oak , _ ii _ Sr es 
Emitter transition capacitance, C,, = Cea/ (Wy Vie) 4. Ro 


and 

Sesto : > _ Te eee 
Collector transition capacitance, Cac = Cao/ V5 Vo) 4.34 
maere n = 1/2 and n = 1/3 for an abrupt junction and a Itnear 


e@eaaed junction, respectively. CEO and Coo are determined in 

the same manner as explained in previous section (Section E). 
Most transistor manuals give the normal and reverse 

common emitter current gains hEE(N) and heecry: The relation- 


Ship between the four gain factors is given by 


Neecny = Gn/ Ol - oy) Oy = Regu)’ * Been?) 
4.35 
MEE(I) ~ %r/ 41 - @y) a7 > Meet)? 2 ey? 
The diffusion capacitances, CHE and Cyc are given by 
Choe = Slee Tp /my, AE See 
and 
Cyc = Olop To /Me re 
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where Tp and Ta are emitter and collector time constants, 


mespecetively. These time constants can be evaluated by the 


Kuno's Method that has been mentioned in Section E. 


= 


CrE Cre 
Cyr Cnc 
ALE B Ri ¢ C 


Ew, ; SB, 
(3 ie 
SS, I a, 
a ad oe * 
7 - 
ee ——— ———_—__—__—-- >} 
Emitter Base Collector 


Figure 4-14, n-p-n Transistor Model for Computer Analysis 


Most o£ the transistor parameters are smppiulcdma, 
he manufacturer. However, if it becomes necessary to deter- 
mine the parameters, Refs. 18, 19, 20 and 21 are very useful 
momeeevaluation. | 


3. Modified Model 


hie temperature GevendenGeme gemma S Pomema Game tomer was 
discussed and some assumptions were stated in Chapter III. Addi- 
tional assumptions are given below: 


(a) ot and a; are approximately the same at temperatures 


ieana (1 + Ai). 
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(b) AT << T, 

(c) TE and tT. are the same at temperatures T anc (T+ AT),. 
(d) Ter = T-clexp(6V_/m,)] and the same approximation 

is for Loe 


(e) Temperature dependence af Rp and Ry is as: explained: 
pare cChapter III. 

Equation 3.17 was written for callector functiom dinde: 

Saturation current, and a similar expression cam be written 


for emitter jucntion diode saturation current as 
TsoT = Iggll + A_C4T)T 


44, 3B: 


Tey * Xe (ATT.1 


EO EO 


where I is emitter junetiom saturation current. at 


EOT 
temperature (T + AT) while ye is at temperature 7... the 
expression for rag Gh given by Eq. 47,51 for a specitied 
memoeerature. If Teor ep Gist Oa eis aegis wesli ible (ess Ic at 


temperature (T + AT), the following expression can be written 


Test = Tgot/( - yey) 


= [lio + dp (AT) Teg lf — Oy) 


lest = Leg/O - Gyoq) * [Ag Tig l/l ~ ayos, 


I 4-. 39 


7 Mag 7 Ma Sec 


ES le 
where Kp = A, (AT). Equation 4.18 was the defining equation 
ffor the factor B. A similar definition can be used for 
emitter junctions as follows: 


Be. = eV./mpk 
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For the derivation of lepr,, we first express the ter 


"exp [B./ (T + AT)]" in @ more convenient form es followss:: 
exp[B,./(T + AT)] = exp(B,/T) exp[-B,(AT)/T(T + AT)] 
= L, exp(8,/T) 
where L, = exp[-B,(AT)/T(T + AT). 
Thus for a change in temperature, AT ,, the following: 


approximate expression for the forward emitter current is: 


obtained: 
ler ~ les [exp (B_/T) ILE 


(Ipc + Ee Ipg) Lexp(B,/T) IL 


Incl 55 K.)L, exp (B./T) 
= Ter (1 + Kp) Ly 


= Me ler 4.40 


where M, = Lp (1 rok 


E BE: 


Another definition given in Section E was Eq. 4.17 
winneh is repeated here for emitter junctions as fellows: 
Ay. = et, /Mpk : 
Then Cyr can be expressed as Cy, = ApI.,/T. Therefore, 


emitter diffusion capacitance at temperature T + AT becomes 


C = A 


pet ~ Ag leer’? * 47) 


Ny. C 


DE 


where Ne = M,T/(T ne 
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As stated before, az large change im h,,.,, however:,, may- 
reflect only a small change in @ amd im many insstancess can 
be neglected in this treatment. Therefore the current. 
generator (a leg) in Fig. 4-14, where the temperature varia- 
tion is not considered, can be writtem a (oy Doe )) with: 
temperature change (AT), so that 

Qylenp = aM. Lex 44, A2- 
in which Eq. 4.40 has been applied. 

One may follow the same derivations given for the: 


collector junction and obtain the following: results: 


*CFY, 7 'q “CF 


Cnocrp = Nc © 


te eee * 
Same manipulation can be done for obtaining 


temperature-caused current generators and capacitors as 


in the diode case and then the following are obtained: 


lepp = Me 


crt > Icr * Icp Me - 1 
DE ee cee DE DE 


= C agen 


pp * Cpg SNe - DD 


Cpcr = Spc * Spe Se - VE « 


Now the modified computer model for transistors can be drawn 


as shown in Fig. 4-15. 


| 2 
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a 
d Ve — es Oe 
Emitter Base 


Figure 4-15. n-p-n Temperature Effected T- 


for Computer Analysis 


G. FIELD-EFFECT TRANSISTOR MODEL 
ieee oeneral 
Although field-effect transistors 
used extensively at present, there 1S no 7 
Mesecribe all such devices. One reason 15 


processes causing various modes of operat: 


another. There are two main types, juncti-: 


SZ 


SLOT 


YQ 











gate FET, sometimes called the "Metal-Oxide Semiconductor" 
(MOS) transistor [Ref. 13]. This would suggest the use of 
at least two different models. A piecewise linear model 
has been developed by Roberts and Harbourt [Ref. 9] that is 
Pesoned for circuit analysis using the ECAP [Ref. 1] pro- 
gram. However, piecewise linear models have some deficiencies 
when applied to general circuit analysis programs because a 
piecewise linear model can be made to approximate the actual 
behavior quite accurately by using several regions. However, 
as more regions are added, more constants must be provided 
and the awkwardness involved in determining the correct region 
fmmoperation. 
fee Models without Température Ettects 

It is possible to make analogy between vacuum tube 
Mele aS a triode) and FET, provided that the cathode is 
replaced by source, the plate replaced by drain and grid 
replaced by gate. Neglecting the bulk resistances and taking 
the barrier capacitances into account, the following equiva- 


fmem@ermeircuit can be obtained. Since the gate junction is 


Gate | Dra ii 





Source Source 


Pueume 4-16. Small-Sicnal JFET Eowivalent €ireuit 


16) 





reverse biased, the Sarees eect resistance cg amd the 
gate-drain resistance Top are extremely igeee amd. ~hearce; 

have not been included in the circuit af Fig. 4¢-I6é.. Deter-- 
Mination of the equivalent circuit parameters will be dis-- 
cussed later. As pointed out before there are two catepories: 
of field-effect transistors, JFET and MOSFET. As the mames: 
moolye a JFEl uses the charatteristics ef a2 reverse—biased. 
mmetion te control the drain-sowrce current, while in. the 
Meee! the gate 18 a métal deposited om an oxide layer ant 

is insulated from the source and drain. Both devices: 

Seerate on the principle of a “elWannel” current comtrolied&k 

by an electronic field. The control mechanisms far the two 
are different, resulting in considerably different characteris-—- 
Eres. The main difference between the two,1is in the gate 
characteristics. The input of the JFET behaves like a 
reverse biased diode while the input of a MOSFET 1s similar 
Beomaesmall capacitor. Based on the above discussion, the 
models shown in Figs. 4-17 and 4-18 are obtained. 

In Fig. 4-17 the source end of the junction_is less 
Strongly reverse-biased or even slightly forward-blased, 
depending on the instantaneous value of Vos: It exhibits- the 
voltage-variable capacitance and conductance properties of 
an ordinary p-n diode. A lumped model of this gate-to- 
source junction can be produced by using the standard built- 
mumodel for a junction diode to simulate these eftects... 


Transconductance and output conductance are simulated by the 
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Gate Dre alate 


\ / : 
Ves Lng ot “ns ps 


¢ @, 


Source Source 


Figure 4-17. Lumped Junction FET Computer Model 


—=— r ae 


Gate | sy. bb eretaisel 
Cn 
Ves Cos 45 DS {js “ns 
Source Source 


Figure 4-18. Lumned MOSFET Computer Model 


Momlinear current source, Ines which can be modelled by 
pointwise linearization as described later. 

AS Stated before, the gate of a MOSFET resembles 
“egeuattcl-plate Capacitor with a thin dielectric layer. 
Oitematelectric results in a very low gate conductance, and 
Since the thickness is constant, the distributed capacitance 
is uniform along the channel and independent of Vcc. This 
Capacitance can be split between the source and drain ends 
of the channel and both lumped capacitors simulated by 


meandard Capacitor model. 


oO 





For more general modelling (including large-signal 
behavior) by pointwise linearization, the functional relation- 
ship is not known in general but can be approximated emipri- 
cally for a given device by measurements, if the desired 


gegree of accuracy so required. 


= 


Solace: Ula) Ae nieleisrosese 
Ree@aO lien liane wen wei ¢ 
Showing Pointwise 
Panera ni gate on 

Procedure 





Transconductance g is the slope of,the transfer 


characteristic curve, Fig. 4-19. At a point (Vic, Ino) the 
tangent line is given by 
Ing = 8m (es ~ Yes) * Ips aes 


where 
Inc in the pinch-off region iS expressed as in Eq. 4.44. 


7 ¥ n 
Ins = Inocll (Vec/Vp) ] 4.44 


Mitere nN 1S a constant for a given unit, whose value usually 
lies between 1.7 and 2.7 (for many units n = 2 is a good 
assumption), and I,., is the drain current when Vcc = 0. 
The determination of the necessary model parameters will be 


given in the following section. 
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3. Determination of Model Parameters 
Capacitance determination: For the model represented 
MmocCeLLON 2, the commen-source short-circuit input, output, 
and reverse capacitances are given in the following 


eepressions, respectively: 


C = C a> oC 


Iss ~ “pg * “Gs 
Coss = Sng * Eps 
Cass = Sng 


Consequently, the model capacitances may be calculated in terms 


of the short-circuit capacitances as 


Chg = rss 
Cos = Cros - Erpss 
Cns = Coss - SRss 


If the capacitances Croo> Cocgs and Cpog are not Supplied by 
the manufacturer the following measurement method can be used 
to determine the unknown values. Three terminal instruments 
are the most useful, since in addition to measuring the 
Capacitance between two terminals, they are equipped to 
eliminate the effects of capacitance between these two termi- 
mans ala a third terminal. This is shown in Fig. 4-20. In 
Fig. 4-20 the presence of Cy and C, does not affect the 
measurement of C\. In some instruments the "low" terminal 

is grounded and the terminal connected to both Cc, and C, 1s 


called the "guard" terminal. Capacitance measurement jigs 


37 





are diagramed in Fig. 4-21. The limitations on this 
measurement are: (1) the gate signal should usually be 
limited to less than 200 mV, (2) frequency should be kept 
between 500 KHz and 5 MHz. Below 500 KHz the chokes must 
be large, and above 5 MHz series resonance and some other 
detrimental effects occur. 

Conductance Measurement: The pinch-off voltage, Vin? 
eie drain saturation current, Inss> the Giain Conductance , 
Zp» and the transconductance g, are four of the parameters 
meguired in the models presented in this section. For the 


pinch-off region, the transconductance is given by Eq. 4.45. 


gn = (-2I)oo/V,) [1 ~ (Vec/Vp) ] CB 
ee Vac = 0 then g, iS maximum, . 
Some Ingg/Vp 
One 
Up 182 Ings/8max 


imesexperimental circuit for determination of Bea Ts given 


toerae., 4-22. 
ihe determination of gy) can be done by a circuit 
aemomown in Fig. 4-23, in which g, (micromho}) = 100 V, (mV) 


when Vy = er aol t 


Measuring : 
Instrument High 


ke taco 1 x 

lEevenie 

sO L-C Low C 

Meter x ; 
Ground a: 2 





Figure 4-20. Three Terminal “‘easurement 
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(b) 
G * 1 
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MG 
6) All capacitors = 0.lyFarads 
me <4 c 
M : All resistors = 1.0 Megahoms 
All RF chokes = 50 millihenry 
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Febeeene 4-21. (a) Croc Measurement, (b) Coc Vieasti cme nuer 


Gel EG Measurement Jigs 
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muotire 4.22.) Circuit for measurement Of oa 
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20 Kohms 


50 
ohms 
enero Mi 
ive farad | V, : (0) 
Z0 Kohms 1 KHz 
1 Volt 
ue =v | RMS 
IWILIS 160 ome) 
10 farad 
ohms 
rine 4-25. Circume 
for measurement of Ep: 
Current Measurement: : 
i is the easiest determined parameter as far as the 


DSS 


Miesmnine Circuit is concerned. The following configuration 


@ameee set for determination of I 


DeSeeye 


=, = 10 volt 


— 


‘Ww 


Blcure 4-24. Circuit for Measurement of Inog 
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4. Modified Model 

As pointed out in Chapter III, the major effects: of 
temperature on a FET are to alter the average channel comnduc-- 
Mivaty, O, and the gate-channel junction barrier potentral, vw. 
Since o decreases with increasing temperature, the effect. of 
this to decrease the drain current as the temperature rises:.. 
On the other hand, || decreases with increasing temperatures. 
causing the magnitude of the total gate-source potential to bes 
reduced, hence increasing the drain current... 

Next, consider the three terminal FET to be operated: in- 
the saturation region with a constant gate-source emf.. There-- 
eore, I = Ino (59). The differential equatron. for the drain. 


DS 
@urrent is given [Ref. 8] by the following relationshin:: 


H) 


dIno/dT (dT n</ap) (dyp/dT) . (31/30) (do/dT) 


gn (dy/dT) + (I) <o/o) (do/dT) 4.46 
Sree IT yo / 9 = ae and Ing 1S proportional to G vat ea en 
temperature T. A reasonable approximation to the mobility- 
eempetTature function [Ref. Z] is 

ES eye 
where n depends on the impurity concentration. Therefore, after 
necessary manipulation the following equation is obtained: 

(1/o) (do/dT) = -n/T 4.47 
If AT << T then dy/dT can be approximated by +#2.2 mV/°C 
_(depending upon whether the transistor is p-channel or 


n-channel, respectively) at low source current level... For many 


purposes n may be approximated to Z (actually 1t ranges fYOM. acme 


ot 





2.7). Taking the above approximations tmto account thes 
Eq. 4.48 can be obtained for draim cunrent change with. 


respect to temperature. 
dI,o/dT = 2.2(mV/°C) g_ - m I,./T 4.. A8 


It 1s easily seen from Eq. 4.48 that the comdittiom for zero: 


temperature coefficient is approximately giver by 
Lic/8 a SS VATE 30) See 44 49° 


When the change in temperature Is very small compared: to: 
the given temperature, T (AT << T),, Eq.. £..48 cam. bee expressed: 
in the following none 

AT) /AT See em oO) a n. Ij</T 
then 

Al nc = 2ikZ or (An) — sy 2) 2 ae i 4.50 
Therefore, one may represent Eq. 4.50 using two current 
generators with current directions opposing to each other. 
aopeyang the Eq. 4.50 to Fig. 4-17 the FEI computer model 
is obtained for a given temperature T and the temperature 
Giemoae AT << [, as shown in Fig. 4-25. If AT is: large, one 
has to make iterative computations by taking proper AT 
increments. 

The necessary model parameters may be determined as 

explained in the previous section. The currents I, and I, 
caused by the temperature can be expressed in terms- of known 


parameters as 


—" eto 
I, Hi Ui e (AT) 4.5 
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1, = (01 (AT) 4.52 


To be dimensionally correct g. Should -be in millimhos, 


temperature T and (AT) in degrees Kelvin. 


Gate Ob arc mel 
Cy 

Vos ¥ 1hs a2 I, a 1, éo Cys Yps 

Source Source 


Ficure 4-25. Lumped JFET Computer Model with Temperature 
Prrect. 


rd 


The enhancement-type MOSFET has a characteristic 
threshold voltage, Vers Imstedad On the pinch-o1t voltage, Vp» 
but the equations are of the same general form: for 0 < Vos 
- Ver < Vos (pinch-off region) Inc = K(Vec - ad Here K is 
a constant and is a function of channel geometry, doping, and 
eeecerive Carrier mobility. In this section an analysis 
is given for the temperature dependence of Vos when the 
current‘1sS maintained constant and saturation conditions 
are present. It is assumed that interface charge Qog 1s 
independent of temperature and that the effects of a finite 
drain resistance can be neglected. Further, it is assumed 
that the only temperature dependent quantities are Voc (DC 


Pace-source voltage), od, (Fermi potential) and U (effective 
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ie ea2ce Mobmilitim~a. | lhe temperature coefficient as a 


mc tLon Of Ves ee CORVicmie ed aan o |e Dy 
dvac/dt = (ip/g)C-1/u,,) (du, / dT) 
- (do, /dT)t[-2(v¢c - Va)/Vpg] + 1} 


PIilecewenewrtrsot anic the second terms are both positive 


3 
it is evident that there may exist a realizable zero tempera- 


mime cOcCtficient point. Supposing that the mobility is pro- 


rd.) and hem tinh 2mOme —Cairy OD Gain s(by 


pentlLonal to T 
taking certain known factors and coefficients into account) 


the following expression for temperature coefficient [Ref. 8]. 


dvgco/dt = (3.0/T) (Vee ~ Ve) 


- [(0.605 - $,)/T]{[2(vee - Vp)/Vpo] -1} 


4.54 


Pmpormted Out by Giralt et al. [Ref. 23],°*a good approxima- 
Semon to the behavior 15 given by 


d]Veco|/dT =h - 2.5|veo - Vpl[mv/°C , Aa 


where Vp 1s the measured threshold voltage obtained from the 
waearont line extrapolation of an I,/g. plot on the vag axis, 
and h is a constant whose value is somewhat dependent on the 
meecorate doping but normally lies between 2 and 6. 

tava ecnenuthe equ 4.5o98cean be written as an 
mde. 4.506, 


Alv = h(AT) - 2.5| Voc - V| (AT) mV : 4.56 


Gs| 


On the other hand when temperature changes as much as AT 


Pane 4.45 becomes 
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pst = 8ml (es * l4Ve51) - Yes] * Ips 


= 8, Vos * 8ml4V%Gsl - 8m Yes * Ips 


Sica glAVcc| : 4.57 
This change AVoc can be obtained either from Eq. 4.54 or 
Eq. 4.56 depending upon the given (or known) information. 
The term (g_|4Vec]) can be defined as temperature resultant 


current generator than temperature dependent computer model 


for MOSFET obtained as shown in Fig. 4-26. 


Gate “Drain 
aie 
acs LDS |AVac| Vv 
mM GS 
Vv S , - 


DS DS 


& 
Source Source 


Figure 4-26, Lumped MOSFET Computer Model with Temperature Effect 
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V. CONCLUSIONS 


fee) GENERAL 

It has been shown that temperature variations have 
various effects on different electramic components.. Resis:- 
tance of a resistor increases or decreases with. temperature. 
depending on the material used. Inductance of am inductor 
may or may not change with temperature depending on its physi-- 
cal dimensions and the material used. Capacitance of @ 
capacitor changes slightly with temperature. However.,, this: 
small change can be extremely important im some situations: 
fe.c., as in resonance circuits). Forward dtrode currents: and: 
diffusion capacitance increase with increasing temperature. 
Giteewettect On a transistor 1S quite complex, and some simpli- 
mae assumptions are listed in Chapter IV. The junction 
forward current and diffusion capacitance increased with 
meneasing temperature. Reverse saturation current also 
miereased. Two types of field-effect transistors were 
examined separately. It was found that with the proper 
operating point, the temperature effect could be eliminated. 
Metis Operating point cannot be attained, the drain current 
changes. Therefore, the effects discussed above should be 
applied to the computer models when one attempts to analyze 
a circuit which is composed of electronic components. 

The models proposed here allow simulation of the behavior: 


of several circuit components. For a reasonably accurate 
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Simulation, the parameters required for the mockedss weree few. 
easily determined. The actual circuit im whicm thee component: 
is to be used is a major factor tm the choice of a: modéd: and: 
the use of these models must be done with. an. understanding. of 
the assumed operating regtons and. characteris:tics: of: the: 
components themselves. The parameters af the moded. must: aiso: 
be chosen to suit the particular application,, becauses even- 
tne best models can only approximate the behavior of: the: 
feral devices. 

Emphasis was placed om empirical determinatipvm of: parameters; . 
however, these parameters may also be determined from: the: 
material characteristics and structural geometry Gia tie idem Trees 
being modelled. 

These models are useful in analysis and design of elec- 
mame Circuits. Also, since the model parameters may be 
eeuecwlated from intormation on the construction of the device, 
Mmemsonould be possible to predict the behavior of a type of 
component before it is fabricated. 

The models are general so that they are compatible. with 


Hees Oot the electronic Circuit analysisS programs in use today. 


EPeeeeoeGGESTIONS FOR FURTHER STUDY 
Pee Modiricattonm of Models with High Frequency Eftects 
High-frequency models for the components are more 
complex than low frequency models,. since additional. reactive 
elements must be included. As a. suggestion for: further study, . 


eRe emaoht consider the effects of high frequency eifects 
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me 


upon these components under normal aqperating: temperatures: 
and temperature as a variable. 
2. FET in Triode Regton 
AS@Ldat sasmeeicu rel si> concermed,, ther developments: 
were done for the saturatiom regitan aperation.. Thee other: 
very important application tis the trirade-repviom operation. . 
In triode region, (0 < Yong < Veo 7 Vp)! thes gio- 1s expressed: 


as 
2 
Ens = Clpgs/Vp) [2Vgg - Vp) - 2gs3 |. Soe 


At an arbitrary point (V Inc) am the: output: characteristic: 


DS” 
curves, the pointwise linear approximation. to: thes output: curve: 


is the straight line i + Ihe- wheres gno- 


ns * Spgs “ps ~ Vos? 


is the slope at (V Ins) Gmyen waieth Ea... 5.12. Then. 1t- 


bre, “4 
should be simple to develop models for triode TevU LOIN) elo ye 
Poemous FEIT*s by taking the température effects into account. 
fee oome Otner Device Models 
Ones uChtealsco develop computer models tor additional 

mevrecs depending wpon the methods stated im this- paper: with 
Memperature or high frequency effects in mind. 

(a) Zener/avalanche diode madel: Since the avalanche 

diode is nearly identical to an ordinary diode for 


operating voltages greater than V, (breakdown voltage), 


B 
it should be possible to modify the Fig. 4-9, diode 
computer model to include the avalanche phenomenon. 


(b) Tunnel diode model: The unique properties: of: tunnel. 


diodes which are most often utilized are its negative 
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resistance and its very rapid switching from Vp 

(peak voltage) to a point in the forward conduction 

Bolo. MoOStmrunine Peadlode Circuits are designed around the 
peak point and valley point. Any tunnel diode model should 
accurately approximate these points as well as modelling 
other important properties. It is relatively easy to 
measure the peak point, valley point and the maximum 

Slope in the negative conductance region. There is 

very little resistance to current flow when reverse bias 
iaeapelred to the junctions; 1.e€., with reverse bias, 

the tunnel diode is nearly a short circuit. 

feyeUnijunction transistor: On the emitter characteristics 
the peak point is defined 
by the expression 

Vo =v + av where v 


D BB 


eee by the In = 


curve. The quantity a, 





called the intrinsic 
standoff ratio. When I. becomes greater than I, (peak point 
current), the process of conductivity modulation occurs in 
the semiconductor material between emitter and base-l, 
causing the negative resistance characteristic. As I, 1s 
increased well beyond I, there is also a slight conductivity 
modulation of the B, - B 


i Z 
s, A model which will adequately describe the static 


region resulting in an increased 


ie 


behavior of a unijunction transistor can be formulated. 
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(3) Silicon controlled rectifiers: Since sxliticom conmtrod. 
rectifier (SCR) is a four-layer, three—-juncttom device with. 
only three leads, it is difficult to describe the internal. 
behavior. The two-transistor analogy provides: a good 
intuitive fee] for the operation, but rt ts duffitcult: to: 
obtain the transistor parameters since base-li and collector=2- 
are not accessible. [t should be possible tm represent: 

an SCR by empirically relating the oie behavioa~ to- 
terminal characteristics even the physical processes: are 
somewhat more complex than for some of the devices: dis-- 


cussed previously. 
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